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Abstract 

We use the membrane paradigm to analyze the horizon dynamics of a uniformly boosted black 
brane in a {d + 2)-dimensional asymptotically Anti-de-Sitter space-time and a Rindler acceleration 
horizon in (d + 2)-dimensional Minkowski space-time. We show that in these cases the horizon 
dynamics is governed by the relativistic CFT hydrodynamics equations. The fluid velocity and 
temperature correspond to the normal to the horizon and to the surface gravity, respectively. 
The second law of thermodynamics for the fluid is mapped into the area increase theorem of 
General Relativity. The analysis is applicable, in general, to perturbations around a stationary 
horizon, when the scale of variations of the macroscopic fields is much larger than the inverse of 
the temperature. We show that the non-relativistic limit of our analysis yields the incompressible 
Navier-Stokes equations. 

PACS numbers: 04.70.-s, 47.10.ad, 11.25.Tq 



I. INTRODUCTION 



The hydrodynamics of relativistic conformal field theories (CFTs) has attracted much 
attention recently, largely in view of the AdS/CFT correspondence between gravitational 
theories on asymptotically Anti-de-Sitter (AdS) spaces and CFTs [l| (for a review see Q]). 
Hydrodynamics is an effective description of the long distance field theory dynamics and 
applies under the condition that the correlation length of the fiuid Icor is much smaller than 
the characteristic scale L of variations of the macroscopic fields. The AdS / CFT correspon- 
dence suggests that the long wavelength dynamics of gravity provides a dual description of 
the CFT hydrodynamics. 

It has been shown in [sl that the {d + l)-dimensional CFT hydrodynamics equations 
are the same as the equations describing the evolution of large scale perturbations of the 
{d + 2)-dimensional black brane. The derivation of this result parallels the conventional 
derivation of hydrodynamics equations from the Boltzmann equation where the "thermal 
equilibrium" solution is the boosted black brane. Thus, the equations of gravity play for 
the hydrodynamic equations of a strongly coupled CFTs the same role as the Boltzmann 
equation plays at a weak coupling. 



xydrodynamics leads to the non- 
61. Since we can obtain the NS 



The limit of non-relativistic macroscopic motions in CFT 
relativistic incompressible Navier-Stokes (NS) equations j^, 
equations in the non-relativistic limit of CFT hydrodynamics, the AdS / CFT correspondence 
implies that these equations have a dual gravitational description, which can found by taking 
the non-relativistic limit of the geometry dual to the relativistic CFT hydrodynamics Q]. 

In this duality picture the dynamics of the fiuid entropy of the field theory at the asynip- 
totic boundary can be expressed directly in terms of the black brane horizon geometry |7|. 
This behavior of the horizon is reminiscent of the membrane paradigm in classical General 
Relativity, where the dynamics of the black hole event horizon is analogous to that of a 
fictitious fiuid j^, y, 10|. Indeed, the membrane paradigm has already been an important 
tool for calculating transport coefficients of the boundary gauge theory [ll, 12, 13, 14, 15 1. 



In 



16j we used the membrane paradigm formalism and an expansion in powers of the 



Knudsen number Icor/L to show that the dynamics of a membrane defined by the event 
horizon of a black brane in asymptotically AdS space-time is described by the incompressible 
Navier-Stokes equations of non-relativistic fiuids. Moreover, the analysis performed in [l^ 
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holds for a non-singular null hypersurface, provided a large scale hydrodynamic limit exists. 



The purpose of this paper is to generalize the analysis of [16| to the relativistic CFT 
hydrodynamics. Our starting point is an equilibrium ((i + 2)-dimensional solution containing 
a timelike Killing vector field and a stationary {d + l)-dimensional causal horizon. This 
solution is associated with a thermal state at uniform temperature. When a hydrodynamic 
limit exists, we can expand the solution in the neighborhood of the causal horizon in powers 
of Icor/L. We consider two specific examples in this paper: a black brane in asymptotically 
AdS and a Rindler acceleration horizon in Minkowski space-time. Assuming there is no 
singularity at the horizon, we show that at lowest orders in Icor/ L the set of Einstein equations 
projected into the horizon surface is equivalent to the {d + l)-dimensional relativistic CFT 
Navier-Stokes equations. Our results imply that the analogy between horizon dynamics and 
hydrodynamics in the membrane paradigm is in fact an identity in certain cases. Since 
the non-relativistic incompressible NS equations arise in the slow motion limit of CFT 
hydrodynamics 5, 6], we will also obtain results of [3] in this limit. 

The paper is organized as follows. In Section 2 we will briefly review some of the basics 
of CFT hydrodynamics that we will need, including the expansion in the Knudsen num- 
ber and the form of the stress-energy tensor in the ideal and dissipative cases. In Section 
3 we will outline the membrane paradigm formalism and consider the geometry and dy- 
namics of null hypersurfaces and of the stretched horizon. We also clarify the relationship 
between the membrane dynamics in general and real hydrodynamics. In Section 4 we will 
apply the membrane paradigm to a uniformly boosted black brane in a (c? + 2) -dimensional 
asymptotically AdS space-time and show that the horizon dynamics is governed by the CFT 
hydrodynamics equations. In Section 5 we will show that our analysis is applicable to cases 
other than the black branes in AdS by considering the example of Rindler space associated 
with accelerated observers in -|- 2 Minkowski space-time. We will show that the Rindler 
horizon dynamics is also governed by the CFT hydrodynamics equations. In Section 6 will 
consider the non-relativistic limit of our analysis and re-derive the non-relativistic hydrody- 



namics results of 



161] . Along the way we will clarify the relation between the surface gravity 



and the fluid pressure. We conclude in Section 7 with a discussion of open problems. In the 
following, unless explicitly stated, we will use the convention G = c= h=kB = ^- 
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II. RELATIVISTIC CFT HYDRODYNAMICS 



Conformal hydrodynamics in + l)-dimensional space-time {d > 2) is described hj d+1 
fields: temperature T{x) and the {d + l)-velocity vector field = 0,...,d, satisfying 

UfiU'^ = — 1. The stress-energy tensor of the CFT obeys 

d^T'^" = 0, = , (1) 

and the equations of relativistic hydrodynamics are determined by the constitutive relation 
expressing T'^'^ in terms of the temperature and the four-velocity field. The constitutive 
relation has the form of a series in the small parameter (Knudsen number) 

Kn = l,or/L < 1 , (2) 

where /cor is the correlation length of the fiuid and L is the scale of variations of the macro- 
scopic fields. Since the only dimensionfuU parameter is the characteristic temperature of the 
fiuid T, one has by dimensional analysis that Icor ~ f;- The constitutive relation reads 

oo 

r^^(x) = 5^rr(x), Tr^{Kn)\ (3) 

where Tl^^{x) is determined by the local values of and T and their derivatives of a finite 
order. Keeping only the first term in the series gives ideal hydrodynamics, while dissipative 
hydrodynamics arises when one keeps the first two terms in the series. 

The ideal hydrodynamics approximation for T^'^ does not contain the spatial derivatives 
of the fields. The / = term in ([3]) gives the stress-energy tensor that reads (up to a 
multiplicative constant) 

= Ku + {d + l)u^uj\ , (4) 

where r]^,, = diag[-, +, +, .., -f-]. 

The dissipative hydrodynamics is obtained by keeping the / = 1 term in the series in 
Eq. (I3l). In the Landau frame 0, 17, 18 1, that fixes the ambiguity in the form of the stress- 
energy tensor under a field redefinition of the temperature and velocity, the stress-energy 
tensor reads (up to a multiplicative constant) 

T^v = T^+' [Vf.u + {d+ l)u^u,] - 2r]a^,, (5) 



where the shear tensor a^,y obeys a^^u" = and is given by 

1 1 
(^f^u = 2 (^M^^ + ^^"m + UuU'^dpUf, + u^uPdpUy) - -^dau"^ [q^^ + UpuJ\ . (6) 

The dissipative hydrodynamics of a CFT is determined by only one kinetic coefficient - 
the shear viscosity rj. The bulk viscosity C vanishes for the CFT, while the absence of the 
particle number conservation and the use of the Landau frame allow one to avoid the use 
of heat conductivity 13]. The dimensional analysis dictates rj = F{X)T'^ where F{X) is a 
function of the dimensionless parameters that characterize the CFT. For strongly coupled 
CFTs described by an AdS gravity dual one gets F = l/vr. 



III. THE MEMBRANE PARADIGM 



The four laws of black hole thermodynamics are global statements derived from the 
Einstein equations restricted to quasi- stationary perturbations near equilibrium. In the late 
1970's and 1980's Damour [8] and later Price, Thorne, and collaborators j^, [l^ developed 
a very general analogy between the local, non-equilibrium physics of any horizon and a 
fluid membrane. In this picture the horizon fluid membrane is governed by the Einstein 
equations, which correspond to fictitious Navier-Stokes equations with universal shear and 
bulk viscosities. In this section we will ffist review Damour's approach to the membrane 
paradigm, which involves the geometry and dynamics of null hypersurfaces. Although this 
formalism is elegant, in some cases it may be more convenient or conceptually useful to work 
instead with Price and Thome's notion of a stretched horizon, a timelike surface located just 
outside the true null horizon. We will briefly discuss how in a particular limit the stretched 
horizon becomes the true horizon and the two approaches yield identical results. 



A. Geometry and dynamics of null hypersurfaces 

We will consider a (c? + 2)-dimensional bulk space-time M with coordinates X^, A = 
0, ...,d+ 1 with a Lorentzian metric qab- Let H he a {d + l)-dimensional null hypersurface 
in this space-time defined by a restriction of the bulk coordinates F{X^) = 0. The normal 
vector to this hypersurface, = g^^dsF, by definition satisfies 

gABi^i"" = g^^'dAFdsF = . (7) 
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This condition implies that for a null hypersurface the normal vector is also a tangent vector. 

We can choose a set of adapted coordinates so that hypersurface in the bulk space-time 
is given by x'^^^ = r = = const., and denote the remaining coordinates on the horizon 
surface as x^, fi = 0, ...,d. In this coordinate system t' = 0, so — > i^. 

The first fundamental form of the horizon is the pullback of the space-time metric to the 
horizon surface. In the adapted coordinate system, 

= 9ABe^e^, (8) 

where represents a horizon basis. It follows from ([7]) that the horizon metric ■y^i, is 
degenerate: 7^jy^^ = 0. We next introduce the auxiliary null "rigging" vector m^, which 
is everywhere transverse to the horizon and normalized such that m^iA = 1- A vector 
satisfying these conditions is 

m^ = m'^ = 1, (9) 
with all other components zero. Using this vector, one can write a completeness relation 

7ab = 9ab - ^A^B - ^BrriA, (10) 

where on the horizon the {d + 2) tensor ^^b reduces to the degenerate induced metric. 

In order to construct the generalized second fundamental form for a null surface, con- 
sider the space-time covariant derivative Va projected into the horizon using the transverse 
projector 11;^ = ^b ~ rn'^f.B and acting on the normal vector l^, 

n^Vc^^. (11) 

Since I^Ib = (I7j), we have in the adapted coordinate system 

^aV/^ = . (12) 

This implies that V^£"^ is tangent to the horizon and can be expanded in the horizon basis 

V/^ = e/e^ . (13) 

The mixed index object 6^*^ acts as a 'Weingarten map" from the horizon tangent space 
onto itself and therefore is the generalized extrinsic curvature of the horizon. Together, 
these first and second fundamental forms provide a description of the embedding of the null 
hypersurface in the bulk space-time. 



Consider Lie transport of 'Jab along the null normal vector £ which is given by the Lie 
derivative C^'Jab- This expression can be split into its trace part (the horizon expansion 6) 
and trace free part (the horizon shear (t^ab)'- 

e = Va^^ , (14) 

<yAB = 7A7f V(C^D) - ejAB/d . (15) 

Here 7^ is the projector onto the d-dimensional spacelike cross-sections of the horizon trans- 
verse to i"^. The components of the Weingarten map are 

7^S0c^ = cr^^'^E + OlE/d (16) 

e^^i^ = K{x)e'' (17) 

BA^'mBTc = (18) 

Note that in these formulas we have used the covector itla, which is a one-form tangent to 
the horizon, just clS £ IS cl tangent basis vector. k(x) is the surface gravity and, as we will see 
below, Qc is a covector whose components can be associated with a horizon "momentum". 
Eqn. ( ITTI) follows just from the null geodesic equation 



'Bt 



Kixy"^ , (19) 



with the surface gravity as the "non- affinity" coefficient. 

We assume that the dynamics of the horizon geometry perturbations are governed by the 
Einstein equations, which are (with a non-zero cosmological constant A) 

Rab - {l/2)RgAB + A9AB = SttT^^* . (20) 

The Ricci tensor contracted with the normal vector and projected transversely into horizon 
Rab^^^c "^^^ t)e expressed solely in terms in terms of the intrinsic horizon metric and 
extrinsic curvature using a generalization of the contracted Gauss-Codazzi equation for a 
null surface ^|. In our adapted coordinates it has the form 

RABi^'e^ = D^Qf - d,Q. (21) 

Since the horizon metric is degenerate, one cannot define a unique connection compatible 
with it. However there is a well-defined rigged covariant derivative operator on the hori- 
zon Z)^, which is defined in terms of the bulk connection projected transversely into the 



horizon 22] . In the adapted coordinate system it has the form 

De^e^ = T\,ej. (22) 

The right hand side of Eqn. (|2T|) can be expressed as the covariant divergence of a horizon 
stress tensor 

T(^H)': = e/ - (5f;e. (23) 

From f[T6l) - f|T8|) . we see that the horizon stress tensor has the general form 

T^Hf, = f^mj'' + nj^ + a(^)^; + ^^tid - + f^), (24) 

where we have used <d = 6 + k. 

Consider first the component of fl2Tl) along £, which is the contraction of the Ricci tensor 
with £^£^. This yields the null geodesic focusing equation 

RABi^i'' = -tv^e + k{x)o - e'^/d - a^^Vf^) . (25) 

Imposing the Einstein equation, there is no contribution from the Ricci scalar and cosmo- 
logical constant terms proportional to the metric due to ([7j) and we have 

- tVf^e + K{x)e - e'^/d - ^fVf^) - mt'^^'H^i^ = o. (26) 

The other d components of the Gauss-Codazzi equation can be obtained by projecting fl2T|) 
orthogonal to i using the projection tensor 7^. When the Einstein equations are imposed, 
the terms proportional to the metric again do not contribute because by construction 7^ and 
£ are orthogonal. To write this equation Damour splits space and time {t, x^) by introducing 
a horizon basis £ = dt + v^di and di. The coordinate t parameterizes a slicing of space-time 
by spatial hypersurfaces and are coordinates on d-dimensional sections of the horizon 
with constant t. As a result the equations take the form 

- C,Q, - en, = -dMx) + D.a^'^^l + ^^d,9 - 87rTX#*£^ef , (27) 

where = Q f-. Together the focusing equation fl26|) and fl27|) describe the dynamics of any 
null hypersurface. 

Although the horizon system is an intrinsically relativistic system, looks just like a 
d-dimensional non-relativistic Navier-Stokes equation 

C^V, + eV, = -d,p + 2r/D,a(^)^' + ^d,9 - T^^^'i^ef, (28) 



where Vi = —Qi/Sir is the membrane's momentum, p = k/Sh the fluid pressure, rj = l/lGvr 
the shear viscosity, = |^ the bulk viscosity, and £^ef T^** an external forcing term. 
Moreover, using the formula for the expansion as the fractional rate of change in the horizon's 
cross-sectional area A, 



the focusing equation ( !26l) can be written like a non-equilibrium entropy balance law 



where the Bekenstein-Hawking entropy S = A/4 and Hawking temperature T = k/2'k. 
Therefore there is an complete analogy between the dynamics of a null hypersurface and the 
dynamics of a non-relativistic fluid. 

The viscous entropy production term due to the shear does appear here as one would 
expect. It is important to note, however, that the general horizon fluid does not actually 
correspond to a real fluid because it possesses an unphysical negative "bulk viscosity". 
Moreover, the second term on the left hand side of the entropy balance law does not appear 
in hydrodynamics and reflects the general non-local and teleological (as opposed to causal) 
character of a globally defined null surface. 

These discrepancies with hydrodynamics arise because the general membrane paradigm 
formalism is valid regardless of the size of the Knudsen number Kn = Icor/L. For example, 
the membrane equations fl26l) and fl27j) describe the dynamics of a black hole with spherical 
topology in an asymptotically flat spacetime (e.g. Schwarzschild). In this case the correlation 
length of a fluid will scale as Icor ~ T~^, where T is the Hawking temperature, while 
T^^ ~ th, where the horizon radius. Since the horizon is compact, L can be no greater 
than ~ rH- Thus the dimensionless Knudsen number Kn in these cases is of order unity 
and hydrodynamics is not an appropriate effective description. 

B. The stretched horizon 

The above membrane paradigm results can also be obtained via the stretched horizon 
formalism (see . Since this approach employs the familiar formalism of a (c? + 1)- 

dimensional timelike surface we can avoid the mathematical complications of dealing with 
a hypersurface whose normal vector is also a tangent vector. A timelike surface is also 



(29) 




(30) 
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physically advantageous as a boundary since a null horizon is an infinite redshift/blueshift 
surface. 

To start, we imagine that the causally complete region of the bulk space-time M outside 
the general horizon discussed in the previous section is foliated by a set of timelike hyper- 
surfaces with spacclikc unit normal vector n^riA = 1- As the previous section we can use 
an adapted coordinate system so that the foliation is given by surfaces of x^~^^ — r — const. 
and the stretched horizon coordinates are x'^. The induced metric on these surfaces is given 

by 

hAB = Qab - UaUb (31) 
and the horizon extrinsic curvature is 

= /i^Vcn^ (32) 

One can also consider a unit timelike vector field U^Ua = — 1 normal to spacelike d sections 
of the timehke surfaces. The induced metric on these horizon cross-sections is 

Sab = fiAB + UaUb- (33) 

The "distance" between a given timelike hypersurface and the true null horizon can be 
parameterized by the affine parameter (or a function thereof, a) along a congruence of 
ingoing null geodesies. For example, in the previous section, the set of null geodesies with 
tangent vector m"^, and affine parameter r. The true horizon is at a = 0. To approximate 
the true horizon one considers a timelike surface a <^ 1 and takes the limit a — > at the 
end of calculations. The combinations an^ and aU^ are fixed in this limit such that 

an^ (34) 
aU"" ^ e^. (35) 

In the horizon limit some kinematic quantities will diverge like inverse powers of a and need 
to be renormalized so they are fixed in the limit. For example, one can show that 

aKfU^Ue ^ k{x) (36) 
as<islKE - + \^^aO (37) 

where k is the horizon surface gravity and cr^^^ and 6 are the true horizon shear and expansion 
respectively. In contrast, 

KfUA ^ ^i, (38) 
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so the horizon's d-momentum with respect to (t, x*) can be obtained without any a renor- 
mahzation. 

The dynamics of the stretched horizon is determined by the usual contracted Gauss- 
Codazzi equations. When the Einstein equation is imposed we have 

DJ^[s) = T^^tte>B (39) 

where is the intrinsic covariant derivative and the stretched horizon stress tensor is 

T^^^^ = K^"" - h^'^'K. (40) 

Note that this same quantity appears in the hterature as a quasi-local energy- 
momentum stress tensor. When evaluated in an asymptotically AdS space-time it is the 
Balasubramanian-Kraus stress-tensor |23i]. After introducing appropriate counterterms to 
(UV) regularize at the AdS boundary, a oo, this stress tensor is equivalent to the stress 
tensor of the CFT. 

In the opposite (IR) limit, at the horizon, it turns out that fHOj) is ill-defined. In this 
limit the non-degenerate induced metric h^j^^ on timelike surfaces becomes the degenerate 
horizon metric 7^1,. Therefore the inverse of the horizon metric does not exist and there is 
no unique canonical way to raise and lower indices. However, one can avoid this problem 
by always working with the mixed index stress tensor 

Tis); = - 5;K. (41) 

The Kronecker delta is well-defined in the horizon limit and Eqns. ( !36l) and ( !37|) imply 
that in the true horizon limit (HTl) agrees with the horizon stress tensor we defined in ( HM 
from the Weingarten map. It can also be shown in the same limit that the timelike Gauss- 
Codazzi equations yield exactly the null focusing equation fl26|) and the Damour-Navier- 
Stokes equation (1271) . 

In the following sections we will apply this general membrane paradigm formalism to two 
examples where a large scale hydrodynamic limit Kn <^ 1 exists, black branes in AdS space- 
time and a Rindler horizon in flat Minkowski. We will show that in these cases the focusing 
and Damour-Navier- Stokes equations are exactly the relativistic Navier-Stokes equations for 
a real fluid. Hence the analogy between null surface dynamics and hydrodynamics is actually 
an identity in these cases. 
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IV. BLACK BRANES IN ASYMPTOTICALLY ADS SPACE-TIME 



We will first apply the membrane paradigm to a uniformly boosted black brane in an 
((i+2)-dimensional asymptotically AdS space-time. This is a solution to the vacuum Einstein 
equations with negative cosmological constant 

Rab + {d+ l)gAB = . (42) 

The bulk metric of this unperturbed, equilibrium solution in Eddington-Finkelstein (EF) 
coordinates is 

ds^ = —2u^dx'^dr H —u^u^dx^dx'^ + r'^rj^^dx^dx'^ , (43) 

where T is the Hawking temperature and = (7,71;*) (7 = (1 — v'^)^^^'^) is the {d + 1)- 
velocity. Note that the fi index on the {d + l)-velocity is raised and lowered by the flat 
metric 77^,^ and its norm with respect to this metric u'^Ufj, = — 1. The black brane horizon is 
located at = vrT and its normal vector is given hj = and 

r = M^. (44) 

We want to consider perturbations of this black brane horizon parameterized by allowing 
u^{x'^), T(x^). The horizon location is rnix^) and x^ are the coordinates on the horizon 
surface. The resulting non-uniform black brane is no longer a solution to the Einstein 
equation. However if the velocity and temperature are slowly varying functions of x^ (long 
wavelength, long time perturbations) we can solve the Einstein equations order by order in 
a derivative expansion and calculate the corrected metric and stress tensor order by order 
in Knudsen number. In what follows we suppose that u^{ex^) and T{ex^) and use e as a 
parameter to keep track of the number of derivatives. 

Using (H3ll we find the induced metric on the horizon at zeroth order in e is 

ds]j = j^^dx^'dx'' = {'KTf{ri^^ + u^u^)dx^dx'' = [TrTyP^^dx^dx'' (45) 

where P^i, is the projection tensor onto the rf-dimensional subspace orthogonal to u'^. Taking 
the Lie derivative along u we get the 0{e) expressions 

e = d^u^ + dVi (46) 
a(f = {-KTf (P%P^9(„M^) - d,u^P,,/d) . (47) 
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We have defined V = u^d^ and ^ = In T. The horizon shear is equivalent to (up to an overall 
factor) the usual fluid shear in hydrodynamics ([6]). Note that f|T5|) implies 

We can also use to calculate the other horizon geometrical quantities. From = 
P^m^'Va^^, we find = —\a^i where = Vu^. Using u'^'V^u'^ = ku^ we find at zeroth 
order the surface gravity is 

k{x) = 2nT{x). (48) 

With the kinematical quantities defined, we now consider the dynamics of horizon per- 
turbations in a derivative expansion in e. The dynamics are described by the generalized 
Gauss-Codazzi equation fl?Il) in vacuum 

5^6/ - = 0. (49) 

The Weingarten map can be calculated directly from the bulk covariant derivative along the 
horizon fi coordinates 

9/ = V.u^ (50) 

We start by calculating the zeroth order (in Knudsen number) part of this horizon stress 
tensor. Derivatives of this stress tensor should yield the equations of ideal hydrodynamics. 
At this lowest order we find 

9/ = -k{x)u^u^. (51) 

Note that this is consistent with the first term (IMl) since = —u^. Expanding out (H9ll 
gives 

5.9/ + rV©/ - r V0A^ - = 0, (52) 
where, following the definition fl53l) . we define the connection as 

^"f.u = ^9"^{df,gu\ + d^g^x - dxg^,u), (53) 

evaluated on the horizon. Plugging in fj48l) and using the formula F^^ja = dx In we find 

at 0{e) 

- VTuy - Td^,u''uu - Ta^ - dV^Tu^ - d^T = 0. (54) 
Projecting along u'^ yields 

d^u^" + dV^ = 0. (55) 
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This equation is equivalent to the vanishing of the horizon expansion, ^ = 0. Note that we 
could have read off this equation from the lowest 0{6) part of the horizon focusing equation 

m- 

Since we are after the relativistic NS equations we will not introduce an explicit split of 
space and time {t, x*) as done previously in the membrane paradigm, which would correspond 
to hydrodynamic equations written in terms of the d-velocity {u'^ = (7,7^*)). Instead, 
we just project transverse to u with the operator to obtain the other set of d equations. 
The result is 

a, + P;9,e = 0. (56) 

Together this set of membrane equations is identical to the equations of relativistic ideal 
CFT hydrodynamics 

d^T^" = (57) 

with traceless perfect fluid stress tensor given in (j4]). The focusing equation is equivalent 
to the contraction along u, u^dyT^^ = 0, while the second equation is equivalent to the 
projection orthogonal, Pa^dyT'^" = 0. 

An important characteristic of ideal hydrodynamics is that the fluid entropy current is 
a conserved quantity 9^Jf = 0. In our case the horizon expansion is the fractional rate of 
change in the horizon cross-sectional area 6 = Cu^^^/^, so our results imply that at 0{e) 
the area is unchanged. Indeed, using the equivalent definition 

9 = -Ld^(^un = 0, (58) 

and the Bekenstein-Hawking entropy formula S* = ^4/4, we find 

= \^u^. (59) 

As expected, and in agreement with the entropy current derived in j?!], the fluid entropy 
density is proportional to the horizon area density (vrT)*^. 

A. Viscous hydrodynamics 

We now consider 0{£) terms in the generalized extrinsic curvature and the structure of 
the membrane equations to 0{e'^). An immediate difficulty is that a priori the location of 
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the horizon th is modified at 0{e). To compute this correction it seems one would have 
to know the full 0{e) metric (the solution to all the Einstein equations up to 0{e'^) 32|). 
However, using horizon geometric variables, we can make the following simple argument that 
the location is unchanged at lowest order. Since the only scalar horizon variable at 0{e) is 
the expansion 9, imposing the ideal equation ^ = means that th = 'n:T{x^) + 0{e'^). 
Next, using the zeroth order metric (H3|) . and (!50!) we obtain 

e,'^ = -2uTu,u^ - \i,u^ + (T^: + -^Pl^B. (60) 

There are several important things to discuss pertaining to this result. First, the ideal 
equations can also be imposed in terms of 0(6^) in the membrane equations, so the last 
term involving Q above effectively will not contribute at this order. Second, the question 
again arises whether the 0[e) corrected metric can affect the remaining three terms in ( 160|1 . 
If they do, can we even proceed without knowing the details of these corrections? It turns 
out 0[e) corrections in the metric can introduce 0[e^^ corrections to the shear tensor, but 
these can only appear at higher order, 0[e^\ in the membrane equations. On the other 
hand, the corrections to the near-horizon metric will contribute and modify the first two 
terms of fl60l) at the order we are considering. In particular, from flM|) we see that at 0[e) 
the non-zero contributions to the horizon stress tensor are 

V)/ = -/t^'^M.M^ + 27rT(mWM^ - - \,u^ + < + -^P^Q - b^^K^^\ (61) 

Here ml^\ and i'^-^^ are the first order parts of the surface gravity, covector, and null 
normal respectively. 

These variables are associated with various ambiguities in the geometrical description 
that need to be fixed. The first ambiguity arises from the fact that the horizon null normal 
vector is not unique; any vector obtained by an overall scaling of the original one is still a 
null normal vector: 

r f{x)F. (62) 
This freedom means that unlike the non-null cases, the generalized horizon extrinsic curva- 
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ture is not unambiguously defined (and can be non-symmetric). Under this scaling 2l| 



^ f ^m^ (63) 

/t' /(/« + !) In/) (64) 

a' ^ fa (65) 

e' fe. (66) 

At zeroth order in e, = (i.e. / = 1) is the natural scaling fixed by the equilibrium 
solution, but at 0{e) there is again an ambiguity. Since there is a zeroth order k (the 
temperature), flMl) implies that the scaling is equivalent to the freedom to define We 
choose a / = 1 + 0(e) such that the 0{e) correction to k is zero; that is 

K = 2nT + 0{e^). (67) 

Thus, fixing this ambiguity determines k at first order and relates the surface gravity to the 
temperature. 

The second ambiguity corresponds to the freedom in the form of the bulk metric evaluated 
at the horizon, i.e. the choice of {d + 1) of the {d + 2) functions Yb in 

Qab Qab + di^AYB)- (68) 

One can also see this from the fact that since Trf = 1 exactly (Q, the covector is m,^ = Qj^r, 
evaluated at the horizon. We make the gauge choice that 

P;m« = l{2nT)-'a^. (69) 



This sets the 0{e) correction to m^, such that the term in (|60|1 is eliminated. 

The third ambiguity is analogous to the ambiguity in the definition of relativistic viscous 
hydrodynamics. This ambiguity corresponds to the choice in the definition of the fields T{x) 
and u^{x). One choice of a hydrodynamic frame is the so called Landau frame, defined by 
the d + 1 conditions imposed on the first order viscous correction to the symmetric fiuid 
stress tensor T^-^-^ 

u.Ti;;'^ = . (70) 

These conditions are equivalent to the statement that is an eigenvector of the full stress 
tensor. Physically the Landau choice is that T[x) and u^{x) are such that in the local rest 
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frame at each point the fluid momentum is zero and energy density can be expressed in 
terms of equihbrium quantities, without dissipative corrections. 

We will flx this ambiguity in the gravitational description by the requirement that 

n^'T^^)'. = 0, (71) 

as the analog of the Landau frame choice at the horizon. Imposing this condition, we flnd 
the requirement that the flrst order null normal obey 

= 0. (72) 

Note that the frame choice (ITTi) implies the constraint m^£^ = 1 also at flrst order, i.e. 

t^.^u, = m^w'^ , (73) 

where we use the zeroth order relations m|f'' = —u^ and = u'^. 

With these gauge and frame conditions (!67|) . (!69l) . and (I7T1) we can proceed to calculate the 
membrane equations using fl52|) . Imposing the lower order ideal hydrodynamics equations 
where they are applicable, we flnd 

+r(2) ve(°)/ - r(2)\,e(°)A'^ = o. (74) 

The last two terms represent corrections to the connection from the flrst order metric 7^^'', 
which is composed of flrst derivatives of and T(x). Since the zeroth order part of the 
extrinsic curvature is 2ttTu'^u'^, these two terms reduce to 

27iT (u,V\n - {l/2)u\''d,-f^^l'^ . (75) 

The horizon metric must satisfy '~)\^}u^ = 0. Therefore, in terms of the null hypersurface 
geometric variables it has the general form (up to an overall factor) 

7?j~<) + ^7i^ (76) 

In f l75l) contributions from traceless shear will be zero identically and in addition we can 
impose the ideal hydrodynamics equation ^ = 0. Therefore we flnd both terms vanish. 
Contracting the remaining terms in (TMl) with u'^ yields the scalar equation 

d^u'' + dV^=^a,,a^^. (77) 
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One can readily show this equation is just the null focusing equation up to 0{e^). The first 
order focusing equation implied = + 0(£:^). Imposing this result means that the u^d^O 
and 6*^ terms in are O(e^) and 0{e^) respectively. Therefore at 0(£:^), 

kO = o^,a^\ (78) 

Using (1551) this equation can also be expressed as an entropy balance law 



Td^{J^) = \^a^,a^- = 2r,a^,a^'' (79) 



where t] = ir'^ -^T'^/IQ is a shear viscosity. This agrees at 0{e'^) with Eqn. B. 27 in [3] and 



the shear viscosity is the Kovtun-Son-Starinets ll| universal viscosity to entropy density 
ratio rj/s = l/An since s = {ttTY/4. Note also that since the right hand side of the above 
equation is positive definite, the fluid entropy always increases. Therefore the second law of 
thermodynamics for the fluid is mapped into the area increase theorem of General Relativity. 
Projecting with we find the remaining d equations are 

a. + P\d,^ = ^P\{d^(yZ - d a>J. (80) 

The set of equations fl77j) and flHOj) are equivalent to the projections along and orthogonal 
to u of the equations of 1st order viscous CFT hydrodynamics d^T^^ = 0, with traceless 
fluid stress tensor given by ([5]). These results show that a relativistic CFT fluid flow is 
encoded in the dynamics of a black brane horizon in AdS. Specifically, as long as a large 
scale hydrodynamic limit exists (so we can expand in derivatives) 

D.T^H); = d,T; = 0, (81) 

at least up to 0{e^). 

i is interesting to compare our results with the holographic fluid-gravity correspondence 
in sl. In this case the black brane solution near the AdS boundary (or equivalently the 
{d+ l)-dimensional boundary stress tensor) is expanded in Knudsen number. The Einstein 
equations at 0{lcor/L) projected into this timelike surface are the equations of ideal hydro- 
dynamics and act as constraints on boundary data. In order to obtain the dual {d + 2)- 
dimensional solution at 0{lcor/L) one can integrate the remaining "dynamical" Einstein 
equations into the bulk radial direction, subject to the condition of a regular event horizon. 
The horizon regularity condition fixes the boundary stress tensor and constraint equations 
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at next order to be those of viscous CFT hydrodynamics, with the particular shear viscosity 
to entropy density ratio l/An. The procedure can be continued to higher orders in Icor/L in 
the same way and it has aheady been used to derive the second order transport coefficients. 
Note, that a boundary stress tensor with a shear viscosity to entropy density ratio which 
is different than l/An corresponds to a bulk background with a naked singularity (see for 



instance 



In contrast, our local analysis at the horizon is only applicable at the lowest orders in the 
expansion in Icor/L. At this level, we show that the details of the expanded bulk solution 
are not required in order to obtain the ideal and viscous hydrodynamics equations. The 
membrane equations and the implicit condition that the horizon is regular imply a boundary 
fluid stress tensor with the particular shear viscosity to entropy density ratio rj/s = l/An. 
This is consistent with having a regular bulk. 

The horizon dynamics is able to capture the leading long wavelength viscous dynamics of 
the finite temperature CFT on the boundary, but clearly cannot capture the all wavelength 
dynamics. In order to go to higher orders in Icor/L using our method would likely require a 
knowledge of the bulk, possibly via integration out from the horizon into the radial direction. 



V. THE RINDLER HORIZON IN MINKOWSKI SPACE-TIME 

To show that our analysis is applicable to cases other than the black branes in AdS, 
we consider the example of Rindler space associated with accelerated observers in d + 2 

Minkowski space-time. The metric is typically written in the form 

d 

ds^ = -t^^^dr^ + di^ + Y^ dx'^dx'i , (82) 

1=1 

where k is a constant surface gravity. This metric can be obtained from the standard 
Minkowski metric via the coordinate transformation x° = ^ sinh(Kr) and x'^^^ = C, cosh(Kr). 
Therefore these Rindler coordinates cover only a "wedge" of the full Minkowski space-time. 
To a uniformly accelerated observer with worldline ^ = const., the surface ^ = is a causal 
boundary that prevents the observer from an access to the entire space-time. 
To employ the membrane paradigm conveniently, we make a coordinate change 

e = r/K (83) 
x'' = Kx\ (84) 
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where r is an affine parameter along ingoing null geodesies. In these coordinates the metric 
f l82|) has the Schwarzschild-like form 

d 

ds^ = -nrdr'^ + {AnrY^dr'^ + '^^ X] ^x^dxi. (85) 

i=l 

When boosted uniformly in the x^ = (r, x*) directions the metric becomes 

ds"^ = -urUf.Uudx^dx" + (4Kr)"Mr^ + n^P^.ydx^dx" , (86) 

with ((i+ l)-velocity defined as in the previous section. Finally, in Eddington-Finkelstein- 
like coordinates x^ = {t, x*) 

ds"^ = —KTu^u^dx^dx'^ — u^dx^dr + K^P^^dx^dx'^ . (87) 

The fiat Rindler metric is a solution to the vacuum Einstein equation with zero cosmo- 
logical constant 

Rab = 0. (88) 

To perturb this horizon we allow for, as in the black brane example, slowly varying (char- 
acteristic scale L ^ fluid velocity u'^(x) and k,{x). k, can be naturally identified with 



a temperature in the following way. Unruh 



25l | showed that accelerated observers feel the 



quantum vacuum to be a thermal state at temperature T = a/27r, where a the observer's 
proper acceleration. This is essentially a local temperature and can be expressed in a Tol- 
man form T = k/2'kx, where x = V~9tt = is the redshift factor. Given this form, 
we define k = 2txT{x) as the location independent temperature of the system. The non- 
uniform Rindler metric is no longer a solution to the Einstein equation, but as before, one 
can solve order by order for the corrections in an expansion in derivatives of x^ (small 
Knudsen number). 

Using (JHTIl (or (l86l) ). the zeroth order induced metric on the horizon (r = 0) is 

7^. = 47r2T2p^„ (89) 

and taking the Lie derivative along u one finds the horizon shear and expansion have the 
same forms (H71) and (H6ll as in the black brane case. Imposing the membrane dynamical 
equations 

D.T^H)", = (90) 
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and following the analysis of the previous section, it is straightforward to find that these 
equations are again equivalent to the equations of relativistic CFT hydrodynamics 



^yT^''' = (91) 

at least up to 0(£:^). The shear viscosity rj is the same in the black brane case. Therefore 
assuming the Rindler horizon has a Bekenstein-Hawking area entropy density, the shear 
viscosity to entropy density ratio in the hydrodynamics equations is again the Kovtun-Son- 
Starinets ratio rj/s = 

In this case though, the result cannot be understood as mirroring the hydrodynamics 
of a field theory fluid living on an asymptotic boundary of space-time. For example, the 
fluid stress tensor cannot be identified with a boundary CFT stress tensor. Essentially we 
have found that this membrane's dynamics can be re-expressed entirely in fluid mechanical 
language even without a notion of holography analogous to AdS/CFT. The fluid system here 
could be interpreted in terms of the near-horizon degrees of freedom as the vacuum thermal 
state (thermal atmosphere of Rindler particles analogous to Hawking radiation) perceived 



by accelerated observers 26|. 



VI. THE NON-RELATIVISTIC LIMIT AND INCOMPRESSIBLE NAVIER- 
STOKES EQUATIONS 

The hydrodynamics of relativistic conformal field theories is intrinsically relativistic as 
is the microscopic dynamics. In particular, the non-relativistic limit of a relativistic con- 
formal hydrodynamics may not be well defined. Nevertheless, the limit of non-relativistic 
macroscopic motions of a CFT hydrodynamics is definable and leads to the non-relativistic 
incompressible NS equations [s], [y]. 

On the gravity side we implement this limit by considering the slow motion regime where 
the (i-velocity f * is a small perturbation. Temporarily restoring c, the horizon coordinates 
(now split into space and time) become (ct, x*) and the non-relativistic slow motion limit 
corresponds to v'^/c <^ 1. In order to keep track of the different terms we impose the 
scaling dt ~ c~^, f * ~ 9j ~ and we consider c oo. The temperature behaves as 
T(x) = To(l -|- c~'^P{x)), where Tq is a constant and P{x) is the fluid pressure. Note this is 
not the only conceivable slow motion, long distance scaling limit. For example, one could 
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imagine and di scaling differently from each other. However, it turns out our particular 
scaling is a natural choice because it is a symmetry of the incompressible NS equations 6|. 

Consider the scaling limit of the relativistic hydrodynamics/membrane dynamics equa- 
tions derived in the previous sections, fl77|) and flHU]) . At lowest order, 0(c^^), the first 
equation (the focusing equation) reduces simply to the fluid incompressibility condition 

div' = 0. (92) 

The second equation is of 0(c~^) at lowest order. Under the scaling the acceleration a^r — »• 
dtVi+v^djVi, while the ^ derivative reduces to just a derivative of the fluid pressure P. On the 
right hand side, the derivative of the shear tensor contributes (l/2)V^fi after imposing the 
incompressibility condition. The second shear times acceleration term does not contribute 
because it is of higher order. Therefore we arrive at the non-relativistic NS equation 

dtVi + v^djVi + d,P = uV^Vi, (93) 

with kinematic viscosity u = {AttTq)'^. 

In [l^ we derived the same incompressible NS equations directly from the membrane 
focusing equation fl2^ and Damour-Navier- Stokes equation f l27|) discussed above in Section 
3 without needing to know beforehand the fully relativistic hydrodynamics equations. Our 
two derivations are completely equivalent. Damour's special horizon adapted coordinate 
system and choice of horizon basis {i, ei = di) gives the generic membrane equations their 
classic non-relativistic appearance, which earlier allowed us to easily show they are the 



incompressible NS equations. The only difference is that in 16| we required a surface 
gravity 

K = 27rTo{l + P - {l/2)v^) (94) 

while here we must have k, = 27rTo(l + P) in the scaling limit. The discrepancy can be 
traced to Damour's parameterization of fluid velocity i = dt + v^di, while we work with 
I = = (7,7f*). The null normals differ by an overall 7 = 1 + (l/2)f^ + ■ ■ ■ factor. As 
we discussed above in Section 4, an overall scaling in the null normal vectors just amounts 
to a difference in gauge, which at lowest order just affects how the surface gravity is to be 
defined. 
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VII. DISCUSSION 



In the paper we used the membrane paradigm and apphed a Knudsen number expansion 
to analyze the l)-dimensional horizon dynamics of a uniformly boosted black brane in a 
{d + 2)-dimensional asymptotically Anti-de-Sitter space-time and to a Rindler acceleration 
horizon in ((i + 2)-dimensional Minkowski space-time. We showed that the horizon dynamics 
is governed by the relativistic CFT hydrodynamics equations. The fluid velocity and tem- 
perature correspond to the normal to the horizon and to the surface gravity, respectively. 
The second law of thermodynamics for the fluid is equivalent to the area increase theorem 
of General Relativity. 

In the Rindler case, unlike the AdS one, there is no holographic screen at the asymptotic 
boundary of space-time. The result that the horizon dynamics is governed by the {d + 1)- 
dimensional CFT hydrodynamics equations may imply that the near-horizon degrees of 
freedom behave ClS db conformal fluid. 

The derivation of the CFT hydrodynamics equations required a knowledge of the horizon 
embedding and employed a local analysis near this horizon. The results apply to a general 
non-singular causal horizon, as long as there is a separation between the characteristic scale 
L of the macroscopic perturbations and some intrinsic microscopic Ic scale given by the 
inverse of the temperature. The non- singularity requirement was used when contracting the 
Einstein equations in order to obtain the membrane equations. The separation of scales was 
required in order to have a small Knudsen number and a valid derivative expansion. The 
separation of scales does not exist in general. This is the reason why, for example, in the 
general Damour-Navier- Stokes equation fl7r|) . the term di6 does not vanish and leads to the 
assignment of an unphysical negative "bulk viscosity" . 

Taking the non-relativistic limit of our analysis, we obtained the non-relativistic Navier- 
Stokes equations as found in the membrane paradigm approach in [l^. In this way we also 
clarified the relation between the surface gravity and the non-relativistic fluid pressure. 

There are various issues to consider next in the membrane paradigm approach to rela- 
tivistic CFT hydrodynamics. Most important is whether the geometrical formulation of the 
relativistic CFT hydrodynamics as a hypersurface dynamics can provide a new insight to 
the nonlinear fluid dynamics. This is of course relevant also to the non-relativistic limit, 
where we obtained the dynamics of non-relativistic incompressible fluid. 
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Another important question is the construction of the higher order derivative terms in the 
hydrodynamics equations. It is possible that these terms are still encoded in the membrane 
equations at higher orders in e. However, it is not yet clear if the analysis can continue to 



be done locally near the horizon 
bulk space-time will be required 



lypersurface. It is likely a more detailed knowledge of the 



27|. 



One can also consider various generalizations of gravity / fluid correspondence in the mem- 
brane paradigm formalism. As a simple example, since black holes in the presence of an 
electromagnetic field act like a charged membrane j^, 10], it should be possible to derive the 
additional hydrodynamic current conservation equation 







(95) 



for gauge fields in black brane backgrounds 
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Another interesting case to consider is non-conformal hydrodynamics. The conformal 



symmetry is bro 
a scalar field 



ien by non-trivial background matter fields in the space-time, for example, 
291]. Therefore, one would need to consider the scalar field equation near the 
horizon in addition to the membrane Einstein equations with non-zero matter stresses. The 
analysis of the membrane equations should be modified at O(e^), producing a positive bulk 
viscosity dependent on the bulk matter field content. 

Finally, the membrane paradigm approach may offer an alternative route to the hydro- 
dynamics equations in generalized gravity theories, where the Einstein field equations are 
corrected by higher curvature terms 30| and the shear viscosity to entropy density ratio of 
the fluid is no longer simply l/47r 31]. 
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